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ABSTRACT 

We study the dynamics of the Magellanic Clouds in a model for the Local Group whose mass is 
constrained using the timing argument/two-body limit of the action principle. The goal is to evaluate 
the role of M31 in generating the high angular momentum orbit of the Clouds, a puzzle that has only 
been exacerbated by the latest HST proper motion measurements. We study the effects of varying 
the total Local Group mass, the relative mass of the Milky Way and M31, the proper motion of M31, 
and the proper motion of the LMC on this problem. Over a large part of this parameter-space we find 
that tides from M31 are insignificant. For a range of LMC proper motions approximately 3er higher 
than the mean and total Local Group mass > 3.5 x 10 12 A/©, M31 can provide a significant torque 
to the LMC orbit. However, if the LMC is bound to the MW, then M31 is found to have negligible 
effect on its motion and the origin of the high angular momentum of the system remains a puzzle. 
Finally, we use the timing argument to calculate the total mass of the MW-LMC system based on 
the assumption that they are encountering each other for the first time, their previous perigalacticon 
being a Hubble time ago, obtaining M M w + M LM c = (8.7 ± 0.8) x 10 n M Q . 

Subject headings: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: interactions - 
Local Group - Magellanic Clouds 



1. INTRODUCTION 

The three-dimensional velocities of the Magellanic 
Clouds, from proper motion measurements using HST 1 s 
Advanced Camera for Surveys (ACS) and a sample of 
background QSOs (Kallivayalil et al. 2006a, b, hereafter 
Kl and K2; see also Piatek et al. 2008, Kallivayalil et 
al. 2009), are ~ 100 km s _1 higher than those used in 
past theoretical modeling of the Magellanic Stream (MS) 
and, for the LMC, now approach the escape velocity of 
the Milky Way (MW). Consequently, as shown by Besla 
et al. (2007), in a ACDM-based model for the MW the 
Clouds, assuming they form a bound system, are likely 
on their first passage. This claim has re-ignited the dis- 
cussion about the origin of the Clouds. Did they form 
in the outer regions of the Local Group or as satellite 
galaxies of the MW? Proposals for the latter, i.e., that 
bind the Clouds to the MW, include the explicit use of 
smaller velocities, for example, by giving the LMC & 
SMC a common halo (Bekki 2008) as well as Modified 
Newtonian Dynamics (MOND) gravity (Wu et al. 2008). 
A recent study by Shattow & Loeb (2009) argues that 
the past orbit of the LMC is naturally confined within 
the virial radius of the MW if a ~ 14% increase in the 
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MW circular velocity (Reid & Brunthaler 2004; Reid et 
al. 2009) is taken in combination with the lower end of 
the allowed proper motion error-space. 

Apart from uncertainties in the MW potential, other 
puzzles remain, such as the high angular momentum of 
the LMC system (e.g. Fich & Tremaine 1991; Gott 
& Thuan 1978). Even with the modest values for 
transverse velocity (~ 100 km s" 1 ) and apogalacti- 
con distance (200 kpc) used in past MS models, a 
mass of Llmc = 2 x 10 10 M q gives the following for 
the orbital angular momentum, £lmCi of the LMC: 
Llmc ~ Mlmc-RlmcVlmo ~4x 10 14 M©kpc km s _1 . 
This is equivalent to the spin angular momentum of 
the Galactic disk, L disk ~ (2/3)M disk i? disk V r ~ 4 x 
lO 14 M kpc km s _1 , with standard values of Mdisk = 
2 x 10 Mq, i?disk = 15 kpc and V = 220 km s" 1 (Sawa 
& Fujimoto 2005). Since the LMC is in a roughly polar 
orbit, these angular momenta make right angles to each 
other. 

This problem has seen many past iterations: Besla et 
al. (2007) show that the LMC may be on a parabolic 
orbit, and hence the 'angular momentum problem' is 
moot. However, in this case an explanation of why the 
LMC is moving so quickly and on such a different or- 
bit from the other satellites is warranted. For example, 
Fich & Tremaine (1991) comment that the maximum 
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line-of-sight velocity of the MS (—410 km s -1 ; Brims 
et al. 2005) is not only higher than that of the Clouds 
themselves (262 km s _1 ) but also higher than any MW 
satellite within 200 kpc. Morphological studies of the 
satellite populations of the MW and M31 arrive at a sim- 
ilar impasse. Van den Bergh (2006) argues that the fact 
that the Clouds are gas-rich dlrr galaxies and yet occupy 
the small Galactocentric distances usually dominated by 
dSph galaxies may be accounted for by assuming that 
they are interlopers that were originally formed in the 
outer reaches of the Local Group. 

Raychaudhury & Lynden-Bell (1989) contended that 
M31 was close enough to the LMC in its early orbital 
history to cause a significant tide, and that this tide was 
oriented such as to generate the high orbital angular mo- 
mentum of the Clouds. It is worth noting here that their 
analysis was done in the context of external tides on the 
MW-M31 system from more distant galaxies. Building 
on this theory Shuter (1992) and Byrd et al. (1994) con- 
sidered the possibility that the Clouds underwent a close 
encounter with M31 having only recently been captured 
by the MW. Given our new velocities it is especially per- 
tinent to not consider the LMC-SMC-MW system in iso- 
lation and to explore whether the Clouds may have been 
subject to external torques before entering the environs 
of the MW. We revisit this classic problem of quantify- 
ing the torque provided by M31 using the framework of 
the timing argument /two-body limit of the least action 
principle, given two new pieces of information: the new 
proper motions for the Clouds and a new transverse ve- 
locity estimate for M31 from satellite velocities (van der 
Marel & Guhathakurta 2008; hereafter VG08). 

Shattow & Loeb (2009) also included the effect of M31 
in their analysis of whether the LMC is bound to the 
MW. We confirm the results in § 4 of their paper but 
take a different approach here. Within the framework of 
the timing argument /two-body limit of the least-action 
principle, could M31 have generated the high trans- 
verse velocities of the Clouds? There are 4 main ef- 
fects on the Clouds' orbits that bear exploring: 1) the 
effect of total Local Group mass, M tot ; 2) the effect 
of / = Mm3i/Mmw! 3) the effect of M31 proper mo- 
tion; and 4) the effect of varying the proper motion 
of the LMC. There is a large range in values for M to t 
because this is a quantity that cannot be directly mea- 
sured without modeling. Various methods give a range 
of 2 x 10 12 - 5.6 x 10 12 M© (Kochanek 1996; Wilkin- 
son & Evans 1999; Klypin et al. 2002; VG08). The 
quantity / is also hard to constrain observationally and 
we look at a range of values from 0.8 - 2.0 (Einasto & 
Lynden-Bell 1982; Evans & Wilkinson 2000; Klypin et 
al. 2002; VG08). Recent VLBA measurements might fa- 
vor / closer to unity (Reid et al. 2009). The proper mo- 
tion of M31 has not been directly measured, but as men- 
tioned above, indirect estimates exist and indicate that 
it is small (VG08). We investigate the full range of val- 
ues from the literature. Finally, the proper motion of the 
LMC is the only observationally well-constrained quan- 
tity in this analysis (Kl) and we test the whole error- 
space using a Monte-Carlo distribution. 

The primary goal is to assess whether the inclusion of 
M31 in the equations of motion of the Clouds can gener- 
ate a significant torque on the L/SMC orbit. In addition, 
we keep track of the magnitude of the relative tidal force 



on the LMC from the MW and M31 to quantify whether 
there are orbits in which the M31 tide is higher. The tide 
is calculated simply via a double radial differentiation of 
the potential at the location of the L/SMC. While the 
locations of the center-of-mass of the L/SMC are allowed 
to move under the influence of the two more massive bod- 
ies, the potential shapes in all simulations are kept fixed. 
Motivated by the findings of this analysis, we apply the 
timing argument to the MW and LMC, assuming that 
the LMC is on its first passage with the previous peri- 
galacticon being a Hubble time, tn, ago, to calculate the 
mass of the MW. 

In most of the subsequent analysis we present results 
only for the LMC and intend our conclusions to be repre- 
sentative of the LMC-SMC system as a whole. It remains 
a possibility that the Clouds are not in a binary system 
(K2; Piatek et al. 2008). We do not explicitly investigate 
this further as it is still possible to find a binary orbit in 
the error space of the SMC for every given LMC orbit. 
Also, a chance three-body interaction (MW-LMC-SMC) 
seems highly unlikely. It has been speculated, for in- 
stance, that the Clouds are members of a small subgroup 
that was captured by the Local Group (Metz et al. 2009), 
or that fell into the MW at late times (D'Onghia & Lake 
2009). Given their mass ratio (~ 10 : 1) the SMC is not 
expected to play a major role in shaping the orbit of the 
LMC. We thus carry out the analysis assuming that the 
orbital path of the Clouds will be dependent on that of 
the LMC, but we do include the perturbative effect of 
the SMC on the LMC given the SMC mean velocity in 
all our calculations (note, however, that this is an orbit 
in which the Clouds are unbound to each other). We 
await smaller proper motion errors for the SMC (see § 0J 
to lift these assumptions. 

In § [5] we recount briefly the work on the relative mo- 
tion of the MW and M31 and use this as the motiva- 
tion for our Local Group model and methods. In § 3 we 
present results including changes to the orbital history of 
the LMC given these new models. We also present the 
results for MW mass. We conclude in § 01 

2. A MODEL FOR THE LOCAL GROUP 

2.1. The Relative Motion of the Milky Way and MSI 

The Local Group is thought to be decoupled from 
the cosmological expansion and gravitationally bound. 
This is supported by the fact that its two major con- 
stituents, the MW and M31, are seen to be approaching 
each other with a radial velocity of ~ 117 km s -1 (Bin- 
ney & Tremaine 1987). The tangential motion of M31 
has not been directly measured, but indirect estimates 
(VG08; Loeb et al. 2005; Peebles et al. 2001) indicate 
that this tangential motion is small, and thus support 
the assumption that the Local Group is bound. 

Some simple dynamical arguments about the the Local 
Group can be made based on known facts. In 1959, Kahn 
& Woltjer proposed the 'timing argument', in which the 
MW and M31 first moved apart due to general cosmo- 
logical expansion, subsequently reversed paths and are 
now falling into each other under their own gravitational 
attraction, their motions being governed by Newtonian 
dynamics. Assuming a zero angular momentum orbit 
(zero tangential motion) and the known current separa- 
tion and radial velocity, the timing argument requires a 
mass m the Local Group > 3 x 1O 12 M to fulfill this 
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trajectory. An alternative explanation would be that 
the MW and M31 are accidentally passing by, but if 
the galaxies have moved at a constant speed for a Hub- 
ble time, their separation would have changed only by 
~ 117 km s _1 /Hq ~ 1.6 Mpc, i.e., less than the distance 
to the next largest galaxy. So this is not a promising 
alternative (Peebles 1993). 

A more sophisticated treatment of the dynamics of the 
Local Group, pioneered by Peebles (1989; 1990; Shaya et 
al. 1995; Peebles et al. 2001), involves using the principle 
of least-action to calculate orbital solutions from incom- 
plete phase-space information by assuming homogeneity 
of the early universe (see also Goldberg & Spergel 2000; 
Goldberg 2001). Note, however, that the two-body limit 
to this solution is equivalent to the timing argument, 
since the equation of relative motion has the same form 
as the cosmological acceleration equation in the case of 
the evolution of a homogeneous, isotropic mass distribu- 
tion (Peebles 1993). 

The timing argument has been investigated and ap- 
plied widely in the literature (Mishra 1985; Raychaud- 
hury & Lynden-Bell 1989; Kroeker & Carlberg 1991; 
Goldberg 2001). Two recent extensions of the formal- 
ism bear further comment: (1) Chernin et al. (2009) in- 
clude the antigravity effect of dark energy in the Kahn- 
Woltjer model and find a larger Local Group mass than 
in traditional methods of 4.5 x 1O 12 M ; (2) VG08 ap- 
plied the timing argument to obtain a total mass for the 
Local Group, after first deducing a transverse velocity 
for M31 from satellite velocities. The methods they use 
for the latter are a statistical analysis of the line-of-sight 
velocities of 17 M31 satellites, a study of the proper mo- 
tions of two satellites M33 and IC 10 (from Brunthalcr 
et al. 2005, 2007), and an analysis of the line-of-sight ve- 
locities of 5 galaxies near the Local Group turn-around 
radius. A full Monte Carlo analysis of all the uncertain- 
ties involved produces a value of Vt an = 41.7 km s _1 for 
the median of the probability distribution of the Galac- 
tocentric tangential velocity of M31, with la confidence 
interval T4an < 56 km s" 1 (thus the radial orbit applied 
by Kahn & Woltjer is allowed in their solution). The 
inferred la confidence interval around the median Lo- 
cal Group mass obtained from application of the timing 
argument is 5.58±g;|| x 10 12 M Q . 

While this mass is consistent with most applications 
of the timing argument, it is on the high end of what is 
predicted from ACDM motivated galaxy models (Klypin 
et al. 2002; Li & White 2008). Using theoretical con- 
straints to narrow the large space of solutions allowed by 
uncertainties in satellite velocities and in halo extent, the 
Klypin et al. models for M31 and the MW have a favored 
total mass of 2.6 x 10 12 M Q (with a maximum upper-limit 
of roughly twice this amount). The reason for the dis- 
crepancy is unclear, the general wisdom being that the 
timing argument provides too simplistic a view of mass- 
accretion. However, some studies (Kroeker & Carlberg 
1991) have investigated its accuracy in a cosmological 
context and found it an adequate approximation to the 
true mass. 

Our aim is not to build a 'true' Local Group but rather 
to investigate maximal M31 models, i.e., models in which 
M31 can hope to generate large torques on the LMC or- 
bit. Thus the timing argument provides a natural frame- 
work in which to carry out these investigations. 



2.2. Methods 

We estimate the effect of M31 by building on the tim- 
ing argument and considering some representative cases 
of the entire paramater-space: we vary the total Local 
Group mass, M to t, from 2.6 — 6 x 10 12 M Q and also con- 
sider two cases for the proper motion of M31 - the mean 
value from VG08 as well as a radial orbit. Table 1 sum- 
marizes the quantities varied, the step-size of the vari- 
ation or error-space explored. We solve the two-body 
problem for the motion of M31 relative to the MW in the 
same MW-centered coordinate system (X, Y, Z) that we 
used for the Clouds (see Kl and van der Marel et al. 2002) 
and then include M31's potential in the equations of mo- 
tion for the Clouds, knowing the distance between the 
Clouds and M31 at every time-step. In solving for the 
radial orbit we assume a distance modulus of 24.47 (Mc- 
Connachie et al. 2005; Ribas et al. 2005), a radial velocity 
(of M31 relative to the MW) of -117 km s" 1 (Binney & 
Tremaine 1987) and a current location for M31 of (RA 
= 10.68°, Dec = 41.27°, J2000.0). For the non-radial 
case we simply adopt the six phase-space parameters sup- 
plied by VG08. Given the adopted initial conditions, we 
are free to choose the total mass (represented here as 
[i = G(Mm3i + Mmw)) so as to specify the time in the 
past at which the MW and M31 were co-located. Fig- 
ure [1] shows the resulting motion of M31 with respect to 
the MW in our Galactocentric frame for a few different 
values of fi/G and a radial orbit. 

As expected, orbits with total mass n/G < 3 x 10 12 M Q 
do not 'turn around' within tn- In Figure[T]we show that 
as /i increases, the turn around time decreases; the dot- 
ted line shows the past orbit for [ijG = 2.6 x 10 12 M Q , 
the solid line shows the case of fi/G = 4.6 x 10 12 M© and 
the dashed line for fi/G = 5.2 x 10 12 M Q . The non-radial 
orbits have the same qualitative behavior with mass and 
we do not plot them here. As in VG08, all quantities are 
found to vary monotonically with the tangential velocity 
of M31: larger values of tangential velocity lead to larger 
values of total mass, the period of the orbit and the peri- 
centric distance. We show the full Hubble time evolution 
of the orbit in Figure [T] simply as a heuristic exercise. 
There is little value in extending our static analysis to 
significantly earlier times than 5 Gyr ago since stellar 
ages imply that the Galactic disk (and presumably the 
halo) had not been fully assembled at z > 2 (Wyse 2007; 
Cox & Loeb 2008). 

Once the mutual separation of M31 and MW is known 
it is easy to introduce a potential term for M31 in the 
equations of motion for the L/SMC: 

d 2 r L _ ^sQtl-ts |) + 
dt 2 </>mw(| r L |) + 0M3i(| i"L - r M31 |)] + J^-, 

(1) 

where yl is the Galactocentric distance of the LMC, rs 
that of the SMC and Ymsi that of M31. <frs is the po- 
tential of the SMC, (f>MW is the potential of the MW and 
0M3i is that of M31. Fl is the dynamical friction on the 
LMC orbit and Ml is the mass of the LMC. There is an 
equivalent equation for the SMC. In this model for the 
Local Group, the MW and M31 are both approximated 
as isothermal spheres with </>mw,M3i (?") = — V 2 In r, 
^o,mw = 220 km s" 1 and Vo,M3i = 250 km s _1 (Loeb 
et al. 2005). To be consistent with the total mass used 
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in our 2-body formulation we introduce a cut-off radius, 
rh, outside which the density drops to zero (see van der 
Marel et al. 2002): 



P 



_ J Vj/4irGr 2 , if r<r h , 
0, if r > r h . 



The enclosed mass is: 



M(r) 



rV 2 /G, 
r h V 2 /G = M t , 



if r < r h , 
if r > r h . 



(2) 



(3) 



The LMC and SMC are represented using Plummer 
models: 

HAr) = GM ilS /[(r - v L . s f + K\^l\ (4) 

with effective radii (Kl, Kg) of 3 and 2 kpc, and masses 
(M L , M s ) of 2 x 1O 1O M and 3 x 10 9 M Q , respectively. 
The final term, J^*-, accounts for dynamical friction on 
the LMC orbit resulting from its motion through the MW 
dark halo, for which we use the Chandrasekhar formula 
(Binney & Tremaine 1987): 



1? — 



{[erf(X)-^| e -* 2 ]vL,s} 



(5) 

where X — wl,s/^o,mw and p{tl,s) is the density of the 
MW halo at the Galactocentric distance of the L/SMC. 
The strength of this force depends on the mass of the 
L/SMC, and since we are integrating backwards in time 
from the current positions and velocities, the sign of the 
dynamical friction term is that of acceleration (see Fig- 
ure 2). Our results do not depend sensitively on the form 
of this drag because as shown in § 3.1 it does not alter 
the L/SMC orbit in such a way as to bring them signifi- 
cantly closer to M31 in the past. Our conclusions don't 
change even if it is ignored. Further, we do not include 
dynamical friction from M31 in our calculations because, 
as discussed in the next section, the LMC's closest ap- 
proach to M31 in the past 5 Gyr is roughly 500 kpc. 
Thus dynamical friction from M31 is negligible. 

3. RESULTS 
3.1. LMC Orbits 

We use a Monte Carlo scheme to randomly draw 20,000 
LMC proper motion components from the errors in Kl. 
From these initial values we allow the LMC orbit to prop- 
agate backward in time in the Local Group model de- 
scribed in § The SMC proper motion is kept fixed at 
its mean value (K2). We keep track of the orbits that 
take the LMC closest to and furthest from M31. 

As laid out in the introduction, we explore the effects 
of 1) M tot ] 2) /; 3) M31 proper motion; and 4) LMC 
proper motion on the past orbital history of the Clouds. 
It was found that as M to t increases, the orbital period 
and apogalacticon of the LMC orbit is smaller, as one 
might expect. However, M to t does not have any direct 
bearing on the relative importance of M31 in the LMC's 
past. The question becomes, given a value of M to t, what 
are the salient changes to the LMC orbit from varying /, 
the LMC proper motion and M31's proper motion? The 
largest change to the orbit comes from varying the LMC 
proper motion, followed by changes in /. The proper 
motion of M31, by contrast, does not dramatically alter 



the picture. It acts much like M to t '■ larger values of 
M31 proper motion require larger values of M to t to fulfill 
the trajectory of the timing argument. The currently ac- 
cepted error-space of M31's proper motion does not alter 
its relative distance to the LMC dramatically. However, 
the proper motion of the LMC does. 

We summarize the results in Figure O The black solid 
line shows the orbit with M to t — 5.6 x 10 12 M Q , / = 1, 
mean LMC proper motion and mean M31 proper motion 
(from VG08). This can be compared to the black dot- 
ted line which is the orbit obtained in a MW-only Local 
Group (with M MW = 2.8 x 1O 12 M ). Note that this 
MW differs from our 'infinite mass' fiducial model in K2 
in that it is modeled with a cut-off radius and thus fixed 
total mass. We realize this is a larger MW mass than 
indicated by either observations or ACDM modeling but 
emphasize that we are trying to understand the relative 
importance of M31 in LMC dynamics for which M to t it- 
self turns out to be unimportant. Varying M to t pushes 
the LMC orbit to larger or shorter period, but the quali- 
tative picture remains. As for the choice of mass distribu- 
tion, we contend that the treatment of M31 as an isother- 
mal sphere is adequate because, as described below, the 
LMC does not come closer than ~ 500 kpc. As for the 
choice of mass distribution for the MW, we showed in 
Figure 8 in Besla et al. (2007) that the trajectories of 
the LMC line up in projection regardless of whether an 
isothermal sphere or NFW (Navarro et al. 1996; 1997) 
profile are used. Thus, replacing the isothermal spheres 
here with NFW models is equivalent to looking at lower 
values of M to t- the LMC would effectively escape sooner 
from the MW in the past and with lower velocity, but at 
the same time, M31 would be further away. 

The dashed pink lines in Figure [2] mark the bound- 
aries of / = 0.8 (minimal M31 model) and / = 2.0 
(maximal M31 model) for fixed M tot . The green lines 
enclose the space of solutions given by LMC proper mo- 
tion: the LMC orbits with closest approach to M31 
have the mean western LMC proper motion, fiyy -, in the 
range of /iw ~ (> 3c), while the furthest orbits are for 
/J-w + (> 3c). These also naturally correspond to un- 
bound and bound orbits to the MW, respectively. The 
northern component of the proper motion, //jy, does not 
affect the trajectory as dramatically (see Figured]). 

Within a period of 5 Gyr, the closest orbits bring the 
LMC within ~ 500 kpc of M31, while the furthest or- 
bits come within ~ 770 kpc. Since the distance to M31 
is relatively large throughout the duration of our calcu- 
lations and tides drop off as r~ 3 it seems unlikely that 
M31 would have a pronounced effect on the orbits of the 
Clouds, and this is the case for most of the parameter- 
space that we search. However, for a range of LMC 
proper motions fiw — (> 3c) the torque provided by 
M31 can be significant. In Figure [3] we show the clos- 
est distance to M31 reached by the LMC as a function of 
the LMC's [iWi^N, and their corresponding error cir- 
cles. The black dots are for M tot = 2.6 x 1O 12 M , 
the green dots for M tot = 3.5 x 10 12 M Q , the blue 
dots for Mtot = 4.6 x 1O 12 M and the pink dots for 
M^t = 5.6 x 10 12 M Q . The M31 tangential velocity is 
that of VG08. There is a clear dependence on fiw but not 
on jLtjy. For the case of —2.3 < fiw < —2.2 mas yr" 1 and 
Mtot > 3.5 x 10 12 M Q , the final relative velocity between 
the LMC and M31 is smaller than the escape velocity 



5 



of M31, but the relative distance is larger than rh,M3i- 
As Mtot increases the value of \iw must also in general 
increase for this to be the case. Varying / broadens the 
distributions in Figure [3] ( left panel) by decreasing (high 
/) or increasing (low /) the distances of closest approach. 

Figure [U shows the past orbits of the Clouds in aitoff 
projection, centered on the MW. The pink triangle, blue 
diamond and red square mark the current positions of 
the LMC, SMC and M31 respectively. Again, the M31 
tangential velocity is that of VG08. We do not show the 
M31 orbit here because, as in the case of the radial orbit, 
the trajectory in aitoff projection is roughly constant. 
For a given M to t the pink solid line shows the LMC orbit 
trajectory that takes it closest to M31, while the pink 
dashed lines show the corresponding furthest orbits. The 
blue dotted line shows the mean SMC orbit for reference. 
As seen in Figure [H M31 does provide a small amount 
of tidal pull on the LMC orbit compared to a MW-only 
Local Group even for the mean velocities. However, M31 
can only alter the orbit of the LMC if the LMC is moving 
roughly 3 — 4a faster than we think, i.e., it is not bound 
to the MW at all. For these cases the tidal influence of 
M31 is larger than that of the MW (solid pink line) and 
the final relative velocity between the LMC and M31 is 
lower than the escape velocity of M31, even though the 
relative distance is larger than r^,M3i- If we do allow 
our calculations of these particular orbits to run for tu 
(not shown here), the LMC eventually becomes bound 
to M31. However, accounting for the hierarchical build- 
up of the Local Group, this analysis becomes inaccurate 
past 5 Gyr. The plot looks qualitatively the same when 
varying M tot . 

Figure [5] quantifies the strength of the tidal force ex- 
erted on the LMC by the MW and M31 as a function of 
time. The tidal force is calculated by twice differentiat- 
ing the gravitational potential of the galaxy in question 
at the position of the LMC, even though strictly the tidal 
field is a traceless tensor that requires tangential deriva- 
tives as well. The radial derivative is proportional to the 
tide only for a perturbed body of constant size (see Gar- 
diner & Noguchi 1996) but is instructive here given our 
methods. Figure \5\ (left) shows the relative tidal forces 
on the LMC orbit that is furthest from M31 (bound to 
the MW) , and on the right for the LMC orbit that comes 
closest to M31 (unbound to the MW). The dotted line 
shows the tidal force exerted by the MW and the dashed 
lines shows the same for M31. The M31 tide takes over 
for orbits moving 3c faster than the mean, but is insignif- 
icant for lower values of LMC velocity. 

3.2. MW Mass 

If the LMC is indeed on its first passage about the MW 
then we can use the timing argument to estimate the to- 
tal mass of the LMC-MW system. This is not strictly a 
correct application of the timing argument because for a 
small galaxy such as the LMC, the details of the neigh- 
boring mass distribution become important and cannot 
be treated as homogeneous and isotropic. Even so, it 
is instructive to see what kind of limits the timing argu- 
ment places on the mass of the MW given the assumption 
that the LMC and MW are encountering each other for 
the first time. Since the current position and velocity 
of the LMC relative to the MW is known, this specifies 
the semi-major axis of the relative orbit via the cubic 



equation 

2n 2 a 3 - n 2 r L a 2 - r L v\ = 0, (6) 

where a is the semi- major axis, n = 2it/T is the mean 
motion, and tl and vl are the position and velocity at 
time zero. If T = t H = 13.7 Gyr (Hinshaw et al. 2008), 
we can solve uniquely for /i = G(Mmw+-^lmc) = n 2 a 3 . 
A Monte-Carlo analysis of the errors in LMC velocity 
gives: 

M MW + Mlmc = (8.7 ± 0.8) x 10 n M Q , (7) 

which is in good agreement with ACDM-based estimates 
for the mass of the MW (e.g., Klypin et al. 2002). It 
is lower than recent observational estimates based on 
VLBA measurements by about a factor of 2 (Reid et 
al. 2009). 

4. CONCLUSIONS 

We have studied the past orbital histories of the Clouds 
in a model for the Local Group whose properties are con- 
strained using the timing argument. This study has been 
motivated by the substantial increase in the tangential 
velocity of the LMC (Kl; Piatek et al. 2008; Kallivay- 
alil et al. 2009) and thus the possibility that the LMC 
is on its first passage about the MW (Besla et al. 2007), 
as well as a new assessment of the tangential velocity of 
M31 (VG08). The timing argument provides a natural 
framework within which to test the effect of 'maximal' 
M31 models on the past orbit of the LMC. The goal is 
to evaluate whether, as in Raychaudhury & Lynden-Bell 
(1989) for instance, tidal torques exerted by M31 are im- 
portant for the LMC. 

We investigate the effects of Mtot, f = Mmzi/Mmw, 
M31 proper motion and LMC proper motion. We find 
that Mtot an d M31 proper motion both act in the same 
way : larger values of M31 tangential velocity lead to 
larger values of M tot , period and pericentric distance of 
the MW-M31 orbit. These in turn affect the period and 
peri/apogalactic distances of the MW-LMC orbit, but 
within the accepted error-space, M31 proper motion does 
not significantly affect the trajectory of the LMC orbit. 
M t ot does affect the trajectory of the orbit if M31 is suffi- 
ciently massive (Mtot > 3.5 x 1O 12 M ) and if the western 
component of the LMC proper motion, \iw is increased 
by roughly 3cr. The amount of this increase depends on 
Mtot- However, for the rest of the explored parameter- 
space, the influence of M31 is negligible. Furthermore, it 
is unlikely that we have underestimated the proper mo- 
tion of the LMC by this amount, and thus conclude that 
the M31 tide is likely insufficient in generating the high 
angular momentum of the LMC orbit. In other words, 
the angular momentum problem remains. 

Finally, motivated by the fact that the LMC could 
be on its first passage about the MW, we calculate the 
implied MW+LMC mass if they were last co-located a 
Hubble time ago. A Monte-Carlo analysis of the errors 
gives: M MW + Af LM c = (8.7 ± 0.8) x 1O U M . Since 
Mmw >> -Mlmc, this quantity is roughly the mass of 
the MW, and is in good agreement with ACDM predic- 
tions (e.g., Klypin et al. 2002). 

It is clear from this analysis that the large tangen- 
tial motion of the LMC cannot be easily explained away 
with manipulations of the Local Group model. Confirm- 
ing this large motion is a high priority. In Kallivayalil 
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et al. (2009) we presented an ongoing analysis of the 
proper motions of the Magellanic Clouds using a third 
epoch of WFPC2 data centered on background quasars. 
The results so far as consistent with those presented in 
Kl. At present the RMS error in the position of the 
quasar is roughly 3 times as large for WFPC2 as for 
ACS. However, with an improved method to deal with 
charge-transfer efficiency and magnitude-related effects, 
and with the increase in time-baseline from 2 to 5 years, 



we expect final error bars for the proper motions that 
are smaller by a factor of ~ 2 from Kl. This, combined 
with our understanding of the properties of the MS, will 
allow us to better constrain the orbit of the Clouds and 
make more specific predictions as to their origin. 

We would like to thank Ed Bertschinger, Paul 
Schechtcr and TJ Cox for useful discussions. 
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TABLE 1 
Parameter-space 



Parameter 


Values 


Step-size 


la error 


Mtot (M Q ) 


2.6-6 X 10 12 


0.1 x 10 12 




/ 


0.8 - 2.0 


0.1 




M31 tangential velocity ( km s 1 ) 


0; 42 




< 56 


LMC proper motion (fJ,N^w) ( mas y r_1 ) 


0.44, -2.03 


Monte Carlo dist. 


0.05, 0.08 


SMC proper motion (hn,IJ>w) ( mas y r_1 ) 


-1.17, -1.16 




Fixed 





Fig. 1. — The relative orbit of MW and M31 is shown in our Galactoccntric X, Y, Z reference frame for a few cases of total mass 
/i/G. The dotted line shows the past orbit for n/G = 2.6 X 10 12 M Q , the dashed line for fi/G = 5.2 X 10 12 M Q and the solid line for 
H/G = 4.6 x 10 12 Af Q . 




j r i i i i 

-5 -4 -3 -2 -1 

time (Gyr) 



Fig. 2. — The Galactocentric distance of the LMC when M31 is included in the Local Group model (/ = 1; black solid line) and when it 
is not (black dotted line). For a given Mtot, the pink lines mark the effect of a maximum M31 model (/ = 2) and a minimum M31 model 
(/ = 0.8). The green lines mark the LMC orbit that comes closest to M31, ~ 3<r increase in fiyy (marked — 3<r), and the furthest orbit, 3cr 
decrease in )xw (marked +3<r). 
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Fig. 3. — The closest approach to M31 as a function of the LMC's /iyy (left) and /xjv (right). The black dots show the case of 
Mtot = 2.6 x 10 12 M Q , the green dots show M to t = 3.5 x 10 12 M Q , the blue dots show M to t = 4.6 x 1O 12 M and the pink dots show 
Mtot = 5.6 x 1O 12 M . 




Fig. 4. — The orbital trajectory of the LMC and SMC in aitoff projection for the past 5 Gyr, centered on the MW. The pink triangle, 
blue diamond and red square mark the current positions of the LMC, SMC and M31 respectively, i.e., at t = 0. The pink solid line shows 
the LMC orbit trajectory that takes it closest to M31 in the past, while the pink dashed lines show the corresponding furthest orbits. The 
blue dotted line shows the mean SMC orbit. 
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Fig. 5. — The tidal force (in arbitrary units) exerted on the LMC by the MW (dotted line) and M31 (dashed line) as a function of time. 
(Left) The forces on the LMC orbit that is furthest from M31, and (right) the forces on the LMC orbit that goes closest to M31. 



